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Abstract 

We define a four-state Potts model ensemble on the square lattice, with the constraints that 
neighboring spins must have different values, and that no plaquette may contain all four 
states. The spin configurations may be mapped into those of a 2-dimensional interface in 
a 2+5 dimensional space. If this interface is in a Gaussian rough phase (as is the case for 
most other models with such a mapping), then the spin correlations are critical and their 
exponents can be related to the stiffness governing the interface fluctuations. Results of 
our Monte Carlo simulations show height fluctuations with an anomalous dependence on 
wavevector, intermediate between the behaviors expected in a rough phase and in a smooth 
phase; we argue that the smooth phase (which would imply long-range spin order) is the 
best interpretation. 
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1 Introduction 



A number of discrete spin models with highly degenerate ground states have critical ground 
state ensembles: these include the triangular Ising antiferromagnet with spin 1/2 (Wannier, 
1950; Nienhuis et al. , 1984) or with general spin (Lipowski et al. , 1995; Zeng and Henley, 
1997), the three-state antiferromagnetic Potts model on the square lattice or equivalently the 
six-vertex model (Lieb and Wu, 1972; Kolafa, 1984; Wang et al. , 1989; Wang et al. , 1990; 
Park and Widom, 1989; Fcrrcira and Sokal, 1995), the three-state Potts antiferromagnet 
on the Kagome lattice (Baxter, 1970), coverings of the square lattice by dimers of one color 
(Fisher and Stephenson, 1963) or of two colors (Raghavan et al. , 1997), and the four-coloring 
of the edges of the square lattice (Read, 1992). 

All known models of this type can be mapped, configuration by configuration, to con- 
figurations of integer-valued 'heights' z(x) (van Beijeren, 1977; Blote and Hilhorst, 1982; 
Kolafa, 1984; Zheng and Sachdev, 1989; Levitov, 1990; Huse and Rutenberg, 1992; Kondev 
and Henley, 1995b; Kondev and Henley, 1996); we will call models with this property 'height 
models'. We can interpret z(x) as parametrizing a surface, which (in most of the above 
cases) turns out to be in its rough phase, described by a gradient-squared free energy at 
long wavelengths. By standard methods, this corresponds to algebraically decaying spin- 
spin correlations. (See e.g. Nelson (1983) and Nienhuis (1987); we review the derivation in 
Sec. 2). The exponents may be calculated from the stiffness constant in the long- wavelength 
Gaussian free energy, and indeed the most accurate way to determine the exponents is by 
directly determining this stiffness constant from the mean-square fluctuations of the Fourier 
transform of z(x), as measured in Monte Carlo simulations (Kondev and Henley, 1995b; 
Raghavan et al. , 1997; Zeng and Henley, 1997). It can be seen that the spirit of this 'height 
model' approach is geometrical: spin fluctuations are visualized as undulations of an inter- 
face model, and critical exponents are obtained through computations of reciprocal lattice 
vectors in a high dimensional lattice. 

The height model we introduce here is a new way of generalizing the three-state Potts 
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antiferromagnet (AFM) on the square lattice to general q (and particularly to g = 4). Recall 
that the three-state Potts AFM is in a ground state if and only if the four spins around every 
plaquette are (ABAB), (ABCB), or a configuration related to these by symmetry. The usual 
g-state Potts AFM with g > 3 allows an additional kind of plaquette configuration in the 
ground state, of type (ABCD), and it turns out that this excludes any height representation, 
which was responsible for the interesting critical properties of the g = 3 case. (Indeed, the 
Potts AFM with g > 3 has exponential rather than power- law decay of correlations (Grest 
and Banavar, 1981).) Therefore, our model excludes the additional kind of plaquette so as to 
keep the height representation. We call such a model a 'constrained Potts antiferromagnet' 
(CPAFM). Our motivation to study this model is that the height space (in which the z vectors 
live) has an especially high dimension (five), and hence might exhibit a new universality class. 

In the remainder of this paper, we present a study of the 4-state CPAFM model on 
the square lattice. First, in Sec. 2, we review the basic ideas of height models (using the 
3-state Potts AFM as a specific example). In particular, we explain how to express a spin 
operator in terms of the height correlation function. In Sec. 3 we carry this out for the more 
complicated case of the 4-state CPAFM, presenting the height mapping of spin configurations 
and the height representation of spin operators. Sec. 4 describes the cluster Monte Carlo 
algorithm needed for studying such models. We present the simulation results, and discuss 
their implications for whether the system is in a critical or a long-range-ordered state, in 
Sec. 5. Finally, Sec. 6 summarizes our results in the context of related models. 

In the present paper, we have considered only the zero temperature behavior. At small 
nonzero temperatures, there are rare plaquettes where the constraint rule is violated. Each 
such point corresponds to a vortex-hke defect (a dislocation) in z(x). In other height models, 
such defects are associated with further interesting exponents (Kondev and Henley, 1996). 
Formulas for these exponents (assuming a rough surface) would be straightforward to derive 
for the CPAFM model as an extension of Sec. 3. 
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2 Example height model: The 3-state Potts AFM 



This section is a pedagogical introduction to some basic concepts and methods which are 
vahd for all "height models" . Our of these, we chose the three-state Potts antiferromagnet 
(AFM) on the square lattice (see Figure 1) as the example model for its obvious analogies 
to the more elaborate 4-state model of Sec. 3. (These analogies also led us to select it for 
numerical comparisons). Furthermore, this model had never been explicitly analyzed using 
the geometrical "height-model" language ^, except in the poorly-known simulation paper of 
Kolafa (1984). 

As already noted, a configuration is a ground state if and only if each plaquette is 
assigned cither two or three distinct spins, with like spins situated on opposite corners of 
the plaquette. Our object is to compute the critical exponent governing the decay of the 
correlation function of various operators. The steps to this are, first, to map the spins to a 
height representation described by a gradient-squared free energy, which allows us to compute 
correlations of the height variables (subsec. 2.1 below, in particular Eq. (2.5)); second, to 
relate the spin operators to the height variables (subsec. 2.2 and 2.3 below). The same steps 
would be used for any other height model (Kondev and Henley, 1995a; Kondev and Henley, 
1995b; Kondev and Henley, 1996; Raghavan et al. , 1997) and in particular for the 4-state 
CPAFM (in Sec. 3, below). ^ 

In view of the application to a multi-dimensional height variable in Sec. 3, the most 

important new concepts here are the three lattices associated with height space, defined in 

Subsec. 2.4. In addition, the concept of 'ideal state' is introduced in Sec. 2.2 which eases 
^The triangular Ising antiferromagnet was so analyzed by Nienhuis et al. (1984). 

^The height map described in subsec. 2.1 also works for the q = Z Potts AFM on a simple cubic lattice. 
That model has long-range order in the ground state ensemble (Banavar et al. , 1980) as would be expected 
from the three-dimensional version of our arguments below. At higher temperatures it almost exhibits 
continuous symmetry behavior (Gottlob and Hasenbusch, 1994; Kolesik and Suzuki, 1995). It is intriguing, 
then, to speculate whether its coarse-grained height variable becomes the angle of an XY-model, as does 
happen in the two-dimensional case we describe in this section. 
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(and geometrizes) the determination of the G vector of a given operator. 

The results of this section were aheady derived by den Nijs et al. (1982), but those 
authors had a different focus: they downplayed the exact mapping to height variables, while 
emphasizing the application of results from renormalization group theory and from exact 
solutions of the eight-vertex model. (For our results, a more elementary theoretical level 
will suffice.) However, once the height-space reciprocal lattice vector G is identified for a 
particular operator, the scaling indices r]{G) are found here (see Sec. 2.3) exactly as in den 
Nijs et al. (1982). 

2.1 Height mapping 

For each ground state configuration, we define a height function z{x.) such that the difference 
in height, Az, in going from one site to a neighboring site is -|-1 if the pair of spins are in 
the same order as in A— >B— >C— >A and is —1 if they are in the reverse order. In this way we 
have a height z(x) uniquely assigned to each site x once the height at any particular site is 
specified (see Figure 1(b)). This z^x.) is well-defined (since the sum of Az around any closed 
loop is zero) for any ground state. 

We can coarse-grain the microscopic 2;(x) to obtain the macroscopic height /i(x). We 
now postulate an effective elastic free energy of the form 

F^J dx^K\Vh{x)\', (2.1) 

where K is a stiffness constant. Writing (2.1) in terms of h{q), the Fourier transform of /i(x), 
gives 



q 

where the Fourier transform is normalized as 

Mx) = ^E^'''^Mq) (2.3) 

^ q 

and N is the number of sites in the lattice. All sums over q run over the first Brillouin zone. 



Hence, the partition function is a product of independent Gaussians and the fluctuations 
in h{q) are given by 

Fourier transforming back to real space shows that the height correlation function is 

(|Mx')-Mx)r)~^logr + C, (2.5) 

where r = |r| = |x' — x| is assumed to be large and C is a constant. Models for which 
(|/i(x') — /i(x)p) diverges as r — > oo (as in (2.5)) are termed "rough." Therefore, models for 
which (2.1) is valid (with finite, nonzero stiffness K) are rough, as is the 3-state Potts AFM. 

2.2 Ideal states 

The first step in relating spin operators to heights is to define a set of reference states, called 
'ideal states', such that the spin, (t(x), on site x is a function of the coarse-grained height at 
X. Ideal states are configurations which (i) consist of a periodic pattern of spins, (ii) have 
zero macroscopic height gradient V/i (and are therefore called 'fiat'), and (iii) have maximal 
entropy for local fiuctuations. In general, "ideal states" make it convenient to identify the 
height-space reciprocal lattice vector (see Subsec. 2.3, below) which is associated with the 
locking operator and with spin operators of interest (see Subsec. 2.4). 

The ideal states for the 3-state Potts AFM are the A B/C state (left part of Figure 1(a)) 
and the others equivalent to it by symmetry; in our nomenclature the labels refer to the 
even and odd sublattices, respectively, and a label 'B/C indicates a disordered site with 
spin having equal probabihty of B or C. A first guess for an ideal state might have been 
the AB state (center part of Figure 1(a)), which is clearly periodic and fiat. However, the 
number of states 'close to' the AB state is less than the number of states close to the typical 
A B/C state. (Here distance between states is measured by the number of spin fiips required 
to convert one into the other.) Hence, the AB state is not an ideal state, but rather a 
transition state between A B/C and C/A B states. 
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All the sites in the completely ordered sublattice of an ideal state (e.g. the even sites in 
an A B/C state) clearly have the same microscopic height z — Zo- The values of z on the 
sites in the disordered lattice are ± 1, with the sign depending on the particular value 
of the spin on that site. Since a disordered site is equally likely to have either of the two 
possible spins, it is clear that the macroscopic height of the ideal state is h — Zo- 

We expect a typical ground state configuration to consist primarily of macroscopic do- 
mains of ideal states, separated by thin walls of transition states. The value of h within a 
domain is simply the value of z on the disordered sublattice of the domain. Therefore, the 
macroscopic height difference, A/i, between two ideal domains can be calculated by finding 
the net height change along a path of neighboring sites connecting a site x in the ordered 
sublattice of the first domain with a site x' in the ordered sublattice of the second domain: 

n— 1 n—1 
A/l = ^ Az((7(Xi), (7(Xi+i)) = A^((7i, (7i+i), (2.6) 

1=0 1=0 

where n is the number of bonds in the chosen path connecting x to x', Xj and Xj+i are 
nearest neighbors for i = to n — 1, (jj = o"(xj) for i = to n, xq = x, and x„ = x'. Note 
that n is even if the same sublattice is ordered in both domains, and n is odd if opposite 
sublattices are ordered in the two domains. It is not difficult to show (see Appendix B) that 
the macroscopic height difference between two domains is uniquely determined (modulo 6) 
by the identities of the ideal states in the two domains. Whether two patches of (say) A B/C 
domain have the same height, or heights which differ by ±6 (or more), is not arbitrary but 
depends on the configurations between the domains. For example, if wc follow a path through 
a succession of domains A B/C, C/A B, C A/B, B/C A, A/B C, and back to A B/C, each 
has a macroscopic height h larger by 1 than the previous domain, so that the second A B/C 
domain has a height that is 6 greater than the first A B/C domain. (Sec Figure 2(a).) Once 
the microscopic height on a particular site is set, the ideal state pattern (and in particular 
the spin value on the ordered sublattice) within every ideal-state domain can be identified 
uniquely from its height h. 
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2.3 Correlation exponents 

Our goal is to compute the correlation function of a given local operator 0(x), which is some 
function of the spins cr{x') for x' in a neighborhood of x. To accomplish this, we must first 
convert this operator from a function of the local spin variables to a function of the local 
height variables. Namely, each site belongs to a domain of some ideal state; the spin variable 
on that site is determined by the ideal state label of that domain, and (as emphasized in 
Subsec. 2.2) every ideal state is labeled unambiguously by its height h. 

Since the periodicity of the spin arrangement in each ideal state divides the lattice into 
'even' and 'odd' sublattices, any local operator 0(x) will have similar periodicity and de- 
pendence on which sublattice x is in: 

0{^)^F^^-\h), (2.7) 

where e(x) is the parity of site x (e(x) = 1 if x is in the even sublattice and e(x) = — 1 if x 
is in the odd sublattice). In addition, since 0(x) depends on h only through its dependence 
on the ideal state label, we have 

F^-\h)^F<-\h + a^), (2.8) 

where a± is any integer multiple of 6. 

In light of (2.7) and (2.8) the function Fq^^ can be expanded as a Fourier series: 

OiK)^F^-\h) = J20pe^''\ (2.9) 

G 

where G ranges over the integer multiples of 27r/6. Prom (2.9) we conclude that the corre- 
lations of 0(x) satisfy 

(O(x)O(x')) = J2 0^"^O^f')*(e'(«'^-^''^')), (2.10) 

G,G' 

where h is the height at x, h' is the height at x', and (...) denotes the average over allowed 
configurations. Note that 
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l^^i{Gh-G'h')s^ ^ j^^iG{h-h') ^i{G-G')h's^ ^ l^^iG{h-h')^^i{G-G')h's^^ {2.11) 

where the second equahty follows because [h — h!) is independent of h! . But (e'*^'^"'^')'*') — 
Sg,g': so that substituting (2.11) into (2.10) gives 

(O(x)O(x')) = J2 OpOf'^^e'""^''-'^'^). (2.12) 

G 

Since {h{q)} are Gaussian distributed, so is /i(x) — /i(x'). Moreover, (/i(x) — /i(x')) is clearly 
zero. Hence, we can use (2.5) to show 

(^QiG{h-h')^ ^ ^-iGH\h-hr) ^ cdG'Dr-''^ as r ^ oo, (2.13) 

where 

m = (2.14) 

and c{\G\) = e'^^^/^ > 0. With (2.13), equation (2.12) becomes 

(O(x)O(x')) « 5^c(|G|)0f )0f' ^V-"- 

G (2.15) 

«/(6(x),6(x'))r-^ 

where 

rj = min s r^c : G = ~g~' ^ ^ or O^^ is nonzero > = fjco, (2-16) 

c = c(G'o), and /(e, e') = 2c Re jOco^Go*}- Note that Gq is the smallest positive wavevec- 
tor where the (height space) Fourier transform is nonzero. Also, /(e, e') is not identically 
zero if and only if at least one of Oq^ and 0^^ is nonzero - and this is guaranteed by (2.16). 
Strictly speaking, equation (2.15) is valid as written only if /(e, e') is nonzero for all four 
possible argument combinations. If this is not true, then for (x, x') where /(e(x), e(x')) — 
we should interpret (2.15) as saying (O(x)O(x')) <^ r''^. 

In any case, we see that an effective elastic free energy of the form (2.1) implies that 
correlations in any (nontrivial) operator decay algebraically - i.e., there is quasi-long-range 
order. Moreover, once the stiffness K is known, the decay exponent rj for correlations in 
any operator can be determined from (2.14) and (2.16). Thus, all correlation exponents are 
known once any of them is determined. 



2.4 Specific operators 

Now we are ready to predict (in terms of the stiffness K) the exponents for exphcit operators 
in the 3-state Potts AFM modeL To do this, one should (for a given site x) make a plot of 
the value of 0(x) as a function of the coarse-grained height h that labels the ideal state, as 
shown in Fig. 2; the Fourier components (as in Eq. 9) can then be read off, along with the 
coefficients O^''^ (which are not usually needed). 

In order to define the operators, it is helpful to represent the spin state on each site by a 
unit vector m(x) in the plane, pointing at angles 0, 27r/3, and 47r/3, for spin states A, B, and 
C, respectively. We must first identify the value of Gq corresponding to a given operator. 

The stog'^ererf magnetization operator m''(x) = e(x)m(x) is the order parameter since it 
distinguishes among the 6 ideal states and also has the maximum height-space period of 6, 
as shown in Figure 2(b). Thus, the minimizing wavevectors in (2.16) are Gs — ±-^. The 

TT 

value r]s — — — then follows directly from (2.14). Wang et al. (1989) found rjs numeri- 
cally from Monte Carlo simulations. The exponents were obtained analytically by Park and 
Widom (1989), using the exact mapping of the 3-state Potts AFM to the 6- vertex model 
and the known exact solution of the latter. They found the free energy cost of 'step' bound- 
ary conditions, forcing Ah = ±2 across the system, to be AF = 27r/6L; via conformal 
invariance this implied rjs — 1/3. Given that this is a height model, one could instead di- 
rectly extract the stiffness K — n/dhy inserting the forced tilt Ah/L into (2.1) (which gives 
F = K{Ah)y2L). 

The ideal states are 'ferrimagnetic' - i.e., they each have a net magnetization M; this 
operator has a height-space period of 3 (see Figure 2(b)). Therefore, Gm — ±27r/3 = 2Gs 
and TjM = 4:r]s = 4/3. The ferromagnetic exponent ijm seems to be considered explicitly only 
by den Nijs et al. (1982), who used a Coulomb-gas approach very similar to ours to reach the 
same conclusions. This exponent is relevant and could easily be measured in simulations, 
but this has apparently never been implemented. 

An ideal state can have either the even or odd sublattice as the ordered one. These 
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possibilities can be distinguished by the sign of the operator P(x) = ^^(^) ^ ^a{x.),a{x'), 

n.n.n. x' 

where the sum is over x' which are second neighbors with x. The height-space period of this 
operator is clearly 2 and, hence, Gsub — ±27r/2. The exponent for this operator is irrelevant 
for the (ground state) 3-state Potts AFM model. 

At this point we may identify three one- dimensional lattices which appear in the height 
representation of the three-state Potts AFM. First, the "repeat lattice" in height space has 
period 6, because (as noted in Sec. 2.2) the spin pattern as a function of h repeats itself with 
that period. Secondly, the "equivalence lattice" has period 1, since h —>■ h+1 is induced by 
a global permutation of Potts spins to give a symmetry-equivalent state. Finally, the "height 
space reciprocal lattice" is dual to the repeat lattice and thus has period 27r/6 for this model; 
it is a general fact about height models that the allowed G values (see Sec. 2.3) belong to 
the reciprocal lattice. These lattices were not defined in prior applications of height models, 
since they are rather trivial in one dimension; they begin to become useful when one wants to 
make a systematic comparison between different height models, and they are quite essential 
in dealing with multi-dimensional height models, as in the next section. 
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3 The 4-state constrained Potts AFM 



A more complicated model with a height mapping is the 4-state constrained Potts antifer- 
romagnetic (CPAFM) model on the square lattice (see Figure 3(a)). Now there are four 
possible spins - A, B, C, and D - and, on top of the antiferromagnetic nearest-neighbor 
constraint, there is the added constraint that no plaquette may contain more than 3 distinct 
spin states. In particular, as in the 3-state Potts AFM model, a configuration is a ground 
state if and only if each plaquette is assigned either two or three distinct spins, with like 
spins situated on opposite corners of the plaquette. 

In this section, we will work through all the stages of analyzing this model via its height 
representation. These are (i) discovering all vector components of the height representation; 
(ii) determining the generic elastic theory; (iii) mapping the ideal-states as a point set 
in height space (associated with this is the "repeat lattice"); (iv) evaluating the height- 
space reciprocal lattice vectors, and finally the exponents, corresponding to the operators of 
interest. 

3.1 Height mapping and elastic theory 

Unlike in the 3-state case, height space for the 4-state CPAFM model is multidimensional. It 
has dimension = 5 and can be viewed as the Cartesian product of a 3-dimensional 'even' 
height subspace and a 2-dimensional 'odd' height subspace. For these we shall sometimes 
use the notation (a, b, c)g = (a, b, c, 0, 0) and (a, 6)q = (0, 0, 0, a, b).) 

The change in height in going from a site x to a neighboring site x' is z(x') — z(x) = 
Az = Az^ (8) Az°, where Az^ = {Azi, Az2, A.zs) and Az*-* = (Azi^Az^) are the even and 
odd height change functions, respectively. Here Az^ = Az^((t(x), cr(x')) is a function only of 
the spins on the two sites and does not depend on which site belongs to which sublattice. On 
the other hand, Az'-* = Az*-*(e(x), cr(x), cr(x')) depends also on the parity e(x) of the initial 
site but does not depend on which spin is on which site, i.e. Az*-'(e, cr, a') — Az'-'(e, a', a). 

The identification of the correct height space is nontrivial: for example, had we simply 
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generalized the q — 3 case, we would have omitted the "odd" components, since the 3- 
state height variable is "even". On the other hand, there is also a danger of including 
redundant height components. By redundant, we mean functions of the spins which are hnear 
combinations of other height variables and of bounded functions of the spins. Appendix A 
gives the explicit construction of z(x) confirming the vahdity of its definition. 

All values of Az can be found from Tabic 1, which directly shows the value of Az^((T, a') 
for 12 of the 24 possible argument combinations and the value of Az^(e, a, a') for 24 of the 
48 possible argument combinations. Values of Az for the remaining argument combinations 
follow from Az^(o-, a') = -Az^{a', a) and Az°(e, a, a') = -Az°(-e, a', a), which must hold 
for Az(x) to be single- valued. Figure 3(b,c) shows an even and an odd component of height 
for the example configuration. 

The primitive vectors for the 'odd' subspace, (l,0)o and (0, 1)q, have been chosen to 
have a 120° angle between them as is standard for a triangular lattice; thus the length of 
{x4,X5)q is {xl + xl — x^x^y^^. ^ The reciprocal spa,ce 'odd' vectors, though, are defined so 
that (height) reciprocal-space (1,0)q has a dot product with height-space (1,0)q of unity, 
and with height-space (0, 1)q of zero; and such that reciprocal space (1, 0)o and (0, l)o are 
60° apart, so (g4, ^5)0 has length (g| + + q^qbY^"^- 

As with the 3-state model, z(x) is well-defined for any allowed spin configuration once 
its value at any particular site is specified, and we can coarse-grain the microscopic z(x) to 
obtain the macroscopic height h(x) . The higher-dimensional analog of the elastic free energy 
in (2.1) is 



A detailed consideration of the symmetries of the square lattice was needed to verify that 

the coefficients in (3.1) are diagonal in the gradient indices. (In a general height model, 

they might depend on four indices, two for the gradients in addition to the two for height 
^The symmetry in the definition of the 'odd' heights would be more evident using an alternative repre- 
sentation using three 'odd' components with the constraint 24 + ^5 + = 0. 




a,/3=l 



(3.1) 
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components.) The symmetries of the Potts spin further restrict the matrix K = {Kap) to 
have the form 



















-\K^ 







2 y 



(3.2) 



where the even and odd height subspace stiffness constants, and K*^ , are as yet un- 
determined. (The 'odd' part of the density in (3.1) is just Vh'-'p; the corresponding 
subblock in (3.2) looks nondiagonal due to our choice of primitive vectors for the "odd" 
subspace.) The generahzation of (2.4) is 



(3.3) 



and, similarly, the analog of equation (2.5) is 



fK- 



TT 



log(r) + 



(3.4) 



where Cap is the a(5 component of a constant matric C. The term 'rough' is used to describe 

any height component, /la, for which (|/?.„(x') — /ia(x)p) oo as r ^ cxo. Hence, the even 
(odd) height components in the 4-state CPAFM model are rough if and only if (-^*^) is 
finite. 



3.2 Ideal states 

The ideal states in the 4-state CPAFM model are like those discussed in Sec. 2.2 for the 3- 
state model, except with three possible spins on each site in the disordered sublattice instead 
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of just two^ (see Figure 3(a), left part); for example, the 'A B/C/D' state with A on the 
even sublattice and disordered B/C/D on the odd sublattice. We shall sometimes use the 
abbreviated label 'X' for the disordered lattice, thus this ideal state is designated 'AX'. Since 
there are two choices for which sublattice is ordered and four choices for which spin appears 
on the ordered sublattice, there are 8 ideal states. 

All of our steps here are closely analogous to those for the 3-state model. All the sites in 
the completely ordered sublattice of an ideal state have the same macroscopic height z = Zq. 
The sites in the disordered lattice have height values Zo + Az, where Az has three possible 
values, corresponding to a step from the ordered spin to neighbors with the three possible 
kinds of spin state on the disordered sublattice. Since these Az vectors add up to zero, and a 
disordered site is equally likely to have any of the three possible spin states, the macroscopic 
height of the ideal state is h = Zq. 

The macroscopic height difference between two ideal domains can be calculated by finding 
the height change along a path of neighboring sites connecting a site x in the ordered 
sublattice of the first domain to a site x' in the ordered sublattice of the second domain. 
The generalization of equation (2.6) is simply 

tt — l n — l 

Ah = ^ Az(e(xi), a{xi), (7(xj+i)) = ^ Az(ej, cr,, ai+i). (3.5) 

1=0 i=0 

Note that now we must keep track of the site's parity e(x) in addition to the sequence of 
spins encountered. Altering this parity leaves Ah^ unchanged while causing Ah'^ to change 
sign. 

The differences (3.5) allow us to map out the ideal-states graph for the 3-state CPAFM 

model, of which each node is the h of the corresponding ideal state. Neighboring ideal states 

(those which differ by the minimum number of spin changes) such as AX and XB arc related 
^This is like the intermediate-temperature state in the unconstrained 4-state Potts AFM (Grest and 
Banavar, 1981); however, the low-temperature behavior of that model may be better described by two states 
disordered on the even sublattice and the other two spin states on the odd sublattice (compare Banavar et 
al. (1980)). 



15 



by a Ah vector which is simply the Az vector for the step AB, as hsted in Table 1, and 
similarly for all other pairs. The ideal-state lattice is 3-coordinated (as is visible in either the 
'even' or 'odd' projection) since, for any ideal state, there are 3 ways to make another ideal 
state with the minimal change. Ideal states related by exchanging even and odd sublattices 
turn out to be related by the shiftvectors in the upper half of Table 2, which are simply 
computed from Table 1; the parity e is that of the first site in the string. 

The 'even' space projection of the ideal-states graph is a 'hydrogen peroxide lattice.' (See 
Fig. 1 of Leu et al. (1969), and references therein. The Bravais lattice of this projection is 
body-centered cubic (bcc) with lattice constant 4. The crystallographic basis is 4 points per 
primitive bcc cell, at [0,0,0] (assumed to represent the ideal state 'AX') and at [0,1,-1], 
[1,0,1], and [—1,-1,0].^ All these points are equivalent by symmetry, and the nearest- 
neighbor bonds make 120° angles. Analogously, the 'odd' space projection is a honeycomb 
lattice with bonds of length 

3.3 Repeat lattice 

Before considering operators as periodic functions of the height, we must understand the 
'repeat lattice'. In general this is defined as the set of possible macroscopic height differences 
between domains of the same ideal state in configurations allowed by the model. (Note the 
distinction from the Bravais lattice of the ideal-states graph.) Conversely, the macroscopic 
height difference Ah between domains of different ideal states can never belong to the repeat 
lattice. (This follows since Ah can never be zero in such a case, which can be checked 
exphcitly for our models.) The repeat lattice of the 3-state Potts AFM is just the set of 
integer multiples of 6. 

It can be checked (see Appendix B) that the repeat lattice of the 4-state CPAFM is a 
^The same lattice is also known as the 'Laves graph of degree three' (Coxeter, 1955) and is the backbone 
of the minimal surface known as the 'gyroid' and found in bicontinuous systems (S. Milner, personal com- 
munication). It is interesting to compare this graph with Fig. 4 of Kondev and Henley (1995b), which is the 
three-dimensional ideal-states graph for a different 4-coloring model. 
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centered hyper-tetragonal lattice in 5-space generated by the vectors ^ 

{af = ^(af±af):ze {1,2,3}}, (3.6) 

where 

af = (8, 0, 0)e af = (0, 8, 0)^ a^ = (0, 0, 8)^ 

(3.7) 

aO = (0,12)o aO = (-12,0)0 = (12, -12)^ . 

The vectors (3.6) span a space of only rank 5, as expected since = 5. To represent the 
repeat lattice it is natural to choose the unit cell with af , af , af , af , and af as edges along 
with the 8-point basis {^lajj^ + A2a+ + Asat : ^i, A^, A^ e {0, 1}}. 

3.4 Correlation exponents 

As in the 3-state Potts AFM, the periodicity of the spin arrangement in each ideal state 
divides the lattice into even and odd sublattices, so that any local operator 0(x) has the 
form 

0(x) = F^W(h). (3.8) 

Again, 0(x) depends on height only in that it depends on the ideal state ensemble, so that 

i^o^^Hh) = F^W(h + a^), (3.9) 

where a-*- is any vector in the repeat lattice. We can therefore expand Fq^'' in a Fourier 
series to obtain 

0(x)=F^W(h) = J]0^-)e'«^ (3.10) 

G 

where the sum ranges over all G in the repeat lattice's reciprocal lattice, which we henceforth 

call the 'height-space reciprocal lattice.' Proceeding as for the 3-state Potts AFM, but using 

(3.4) in place of (2.5), shows that 

^These happen to be twice the vectors in the lower half of Table 2. 
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(O(x)O(x')) ^ F(e(x), e(x'))r-'' as r ^ oo, 



(3.11) 



where 




and G{rj) = {G : r^c = v}- In each 



(3.12) 



(3.13) 



GgG(»?) 



case G is in the repeat lattice's reciprocal lattice. Again, equation (3.12) guarantees that 
F(e, e') is not identically zero, and equation (3.11) should be read as (O(x)O(x')) <^ if 
it happens that F(e(x), e(x')) = 0. 

We see from the above and equation (3.2) that the 4-state CPAFM exhibits quasi-long- 
range order if and only if either or is finite. If they are both finite, the algebraic 
decay exponent for the correlations in any local operator can be determined once and 
are known. If either one is infinite, some spin-spin correlation does not decay to zero as 
r — > oo, which is equivalent to the existence of long-range order. 

3.5 Specific operators 

Now we are ready to predict (in terms of and K'^) the exponents rju and 773 for the 
magnetization and staggered magnetization operators in the 4-state CPAFM. We shall rep- 
resent the spins A, B, C, and D as unit vectors m(x) pointing toward the corners of a regular 
tetrahedron. The staggered magnetization m*(x) is defined just as in Sec. 2.4. Determin- 
ing the wavevectors associated with the operators m(x) and m'*(x) is trickier for the 4-state 
CPAFM than for the 3-state Potts AFM because of the higher dimensionality of height space. 
In addition, the optimal wavevector G, the one giving the minimal value of r/c, depends on 
the ratio K^/K^. 

From equations (3.2) and (3.13) it is clear that the optimal G may be purely even (i.e., 
with no component in the odd height space), purely odd, or 'mixed' (having both even and 
odd components), depending on whether / K*^ is large, small, or intermediate in value. 
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The minimal length reciprocal lattice vector G in the purely even, purely odd, or mixed case 
is 



G 



respectively. Each of these is symmetry-related to other Gs with the same length. As it 
turns out, the staggered magnetization operator m*(x) has nonzero Fourier components at 
Qmix g^j^^ qe^ while the magnetization m(x) has nonzero Fourier components at G™'^ and 
GO. 

Eq. (3.13) gives the corresponding exponents, 

(Since rf" and r^^ are independent of the way we parameterized height space, we shall use 
them as surrogates for and in subsequent expressions.) The exponent in (3.12) is 
T]rain = niin [i]^, 1]'^, T]"""^}. lu particular, 

if r]^/r]^ < 1/2, then r?mm = and Vm > Vs = Vmm = ff, (3.16a) 
if 1/2 < rf/r]'' < 3/2, then r]^^ = tT"^ and r]M = m = Vmin = ry™", (3.16b) 
and if 77^/77° > 3/2, then 77^;^ = 77° and 77s > 77M = 77min = ?7°- (3.16c) 

3.6 Locking exponent 

Roughly speaking, there are more states with coarse-grained h near a node than near an 
interstitial point of the ideal-states graph. In effect, an operator /lock(h) is added to the 
effective free energy, a periodic function of h which is negative near each node of the graph. 
This (small) term favors 'locking' of the entire system into a 'smooth' state, which is close to 
one of the ideal states. It is known that locking occurs whenever the corresponding exponent 
77iock is less than 4 (Jose et al. , 1977). 
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This forces our attention to the Bravais lattice of the ideal-states graph. Usually this is 
denser than the repeat lattice, since it includes translations which make the nodes superpose, 
even if they represent different ideal states. Correspondingly, Giockj the smallest reciprocal 
lattice vector of the operator /lock, is typically larger than Gmin- Thus, in the 4-state CPAFM 
case, the 'even' projection of the ideal-states graph has a bcc Bravais lattice, which has twice 
as many points as the simple cubic repeat lattice, and correspondingly has G^^,;^ larger than 
G^ defined in Sec. 3.4. However, when the graph is considered in 5 dimensions, it turns 
out that its Bravais lattice is the same as the repeat lattice, since the 'odd' components 
deviate in the opposite directions from the body-corner and body-center points of the 'even' 
projection. Although the corresponding reciprocal lattice contains G^, there is an extinction 
in the Fourier transform of the ideal-states graph and thus /lock(h) does not have a component 
at wavevector G^. 

We find that the optimal reciprocal lattice vector is 



(Each of these is symmetry- related to other Gs with the same length). Consequently, in 
terms of rf" and 77*^ as used in the preceding section. 




or 




(3.17) 



if Ty'^/r^o < 1/6, then Tyiock = vL^ = 277"^, 



(3.18a) 



if 1/6 < 77''/^° < 15/2, then Tyiock = Co 




(3.18b) 



and if ri^'/rj^ > 15/2, then Tyiock = Ck = ^^r]^- 



(3.18c) 



20 



4 Monte Carlo simulation 



To obtain more information on the 4-state CPAFM we performed a Monte Carlo simulation 
of the model. By averaging appropriate quantities over the ensemble of allowed states, we 
can independently determine the stiffnesses, and K'^, and any correlation exponents. 
The results can then be examined in light of Section 3. 

Although we are only particularly interested in the case q — A,we developed an algorithm 
for simulating the general g-state CPAFM. This permitted us to perform a partial check of 
our code by simulating the q — 3 Potts AFM and comparing our Monte Carlo results with 
the known exact results for that model. 

4.1 Monte Carlo algorithm 

We chose to use periodic boundary conditions in the simulation. Given this choice, it is 
natural to require that the number of rows and columns in the lattice both be even numbers; 
otherwise, there can be no well-defined even and odd sublattices and, hence, the odd heights 
are not well-defined. In addition, there would be no allowed configuration corresponding to 
a single ideal state. 

4.1.1 Initial state 

We performed runs with each of three types of initial state. The simplest choice is one of 
the 2q ideal states. All the spins in one sublattice are set to the same value, while the spins 
in the other, disordered sublattice are chosen randomly from the remaining q — 1 possible 
values. The most important property of this state is that it is 'flat' in height space (the 
coarse-grained height is constant throughout the system). As the system evolves into a 
rough state, the Fourier spectrum (|h(q)p) increases from zero to the equilibrium value, for 
most q. 

A second, more elaborate method is to divide the lattice into an array of smaller blocks; 
each block is assigned a random ideal state (to be filled in as described above), with certain 
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restrictions to ensure that only allowed configurations are obtained for the full lattice. This 
initial state is also macroscopically fiat. 

A third, contrasting kind of initial state for a height model is not macroscopically fiat; it 
has, locally, the maximum possible height-space slope (Kondev and Henley, 1995a; Kondev 
and Henley, 1995b; Raghavan et al. , 1997). It is called a 'roof configuration since it is 
usually composed of two domains of opposite slope. (The system-average slope must be 
zero in order for the states of lowest free energy to be accessible by local updates.) The 
'roof configuration shown in Figure 4 includes 4 domains of periodic states. In each domain 
(considered as a state in the q = 4 model), the odd heights are flat, while the even heights 
vary as rapidly as allowed. In this case, many Fourier amplitudes are large at first and 
approach their equilibrium values from above. Thus, by iterating until the expectation value 
is independent of initial configuration, we can ensure that equilibrium is reached. 

4.1.2 Update move 

An update move must satisfy three conditions: 

(i) it must change any allowed configuration into another allowed configuration; 

(ii) it must be ergodic - i.e., through a sequence of update moves, any allowed configuration 
must be accessible from any other one; and 

(iii) it must satisfy detailed balance. 

In the CPAFM model the first condition is nontrivial: a single spin-flip is disallowed on 

many sites because of the nearest-neighbor and plaquette constraints, and the remaining 

allowed single-spin flips are too few to be ergodic. Thus a cluster-update move is a necessity 

for the CPAFM ground-state ensemble.''' Clearly, within a cluster, the change should be 

a global symmetry of the model, i.e. an exchange of Potts states X and Y. (An allowed 

single-spin flip is a special case of this.) 

''In the 3-state Potts AFM, it is an optional improvement to obtain acceleration (Wang et al. , 1990). 



22 



The update move we used is to 

(1) Pick a site x at random; say cr(x) = X, and randomly choose Y X) to be any of 
the q — 1 values of spin. 

(2) Determine the minimal domain containing site x in which it is allowable to exchange 
X on all sites (we call this the 'XF-cluster' of site x), and do so. 

The geometrical rules for an XF-cluster are simple. We define a hnk (in the percolation 
sense) between sites x and x' if either (i) they are nearest neighbors and have spin values 
X and Y; or (ii) they are second-neighbors, both with spin X or both with spin Y, and 
the other two spins on the plaquette are different. Then the XF-cluster consists of all sites 
connected to x by such links; in this fashion all X and Y spins are divided into disjoint 
clusters, as in a percolation model. Clearly, performing X ■(-^ Y on one endpoint of a link 
and not the other creates a disallowed state, so we must at least update the XF-cluster; it 
is easy to see that it is also sufficient to do so. 

Since the exchange X ■(-^ Y does not change the pattern of XF-clusters, the update move 
is reversible: it can change configuration $ directly into configuration ^ if and only if it can 
change ^ directly into $. Furthermore, the probability of flipping a given XF-cluster in an 
update step is just l/(g — 1) times the fraction of all lattice sites which belong to that cluster 
(since we would have flipped the same cluster if we had initially hit on any of its sites and 
chosen the appropriate second spin value). Hence the rate of $ — given configuration $, 
is the same as the rate ^ ^ given ^; for our ensemble, in which all allowed configurations 
have equal weight, this is the detailed balance property. The only thing left is to verify the 
ergodicity of the update move, which is done in Appendix C. 

In speaking of ergodicity, we should make clear that, in the presence of periodic boundary 
conditions, any local update move conserves the net height difference (tilt) across the system 
and hence is non-ergodic in a trivial sense. However, our cluster update move is nonlocal 
and can change the net tilt. 
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It is interesting to note that, in terms of the height variables z(x), the update move 
simply inverts the even height components on the sites within the updated cluster (the 
cluster boundary is a loop of constant height, as studied by Kondev and Henley (1995a), 
so the height is unchanged along it). This reflection is closely related to the 'valleys-to- 
mountains reflection' algorithm (Evertz et al. , 1991; Hasenbusch et al. , 1992). On the 
other hand, the odd height components on the cluster sites are simply shifted by a constant 
vector. 

4.2 Measurements 

Simulation runs were performed for lattices of various sizes ranging from 16 x 16 to 200 x 200. 
As mentioned above, since the notion of even and odd sublattices is essential in these models, 
the number of sites in a row or column of the lattice was always chosen to be even. In runs 
where Fourier transforms were taken, fast Fourier transform (FFT) algorithms were used 
and, hence, it was also convenient in these runs to have the number of sites in a row or 
column be a power of two. 

Since successive states are highly correlated, any measurements were taken once per 
sampling interval; the duration of a sampling interval (as a number of cluster hits) was 
varied depending on the lattice size. Furthermore, in order to estimate uncertainties in 
the calculated averages, we computed averages and variances using groups of 100 sampling 
intervals. Averages over successive groups were treated as uncorrelated for the purpose of 
estimating uncertainties. 

We made the following measurements. At each sampling time, the height configuration 
z(x) was generated from the spin configuration. Then h(q), the Fourier transform of the 
height configuration z(x), was computed (via FFT). The average of |h(q)p was accumulated 
for every wavevector q. 

We must note that our procedure in generating z(x) is not quite safe. The FFT requires 
that z(x) obey periodic boundary conditions; but if the configuration has a net tilt in height 
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space, the algorithm leaves a 'seam' along the edges of the lattice, across which z(x) has a 
sudden jump. (Recall that the tilt can be changed by the algorithm only in the rare event 
than an update XF-cluster winds nontrivially around the periodic boundaries.) If 'tilted' 
configurations are frequent, this causes a spurious contribution to (|h(q)p) for (only) q lying 
on the X and y axes of Fourier space.^ That would be visible in data plots such as Figure 6: 
those data points would form a separate curve, shifted above the points from elsewhere in 
Fourier space. In fact, the plotted points appear to fall along a single curve. Therefore we 
believe that configurations with height-space tilt are rare and do not affect our results. (If 
the system is in a 'smooth' phase, as discussed in Sec. 5, suppression of tilt is just what one 
would expect.) 

For the 4-state model, additional measurements were made. We computed the total 
magnetization M = ^ m(x) and staggered magnetization M'' = ^ m'^(x), and accumulated 
averages of |Mp and |M*p. Also, we measured the size distribution of the XF-clusters. 



^Such behavior was very obvious in (unpublished) test runs of Kondcv and Henley (1995b). Their 
published data were from simulations that disallowed any update moves that would have changed the tilt. 
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5 Simulation results and discussion 

We will now present and discuss the Monte Carlo data. Our chief interest is whether the 
heights are in a smooth or a rough phase, corresponding respectively to long-range order 
or power-law correlations of the underlying Potts spins. Since the 'even' and 'odd' height 
components are not related by symmetry, it is also possible that one kind could be smooth 
and the other kind could be rough, as occurs in the dimer-loop model (Raghavan et al. , 
1997). 

The 'smooth' scenario is favored by measurements of the height fluctuations, which con- 
tain the most information (Sees. 5.1 and 5.2). On the other hand, direct measurements 
(Sec. 5.3) of the order parameters suggest a 'rough' phase. The evidence will be judged in 
Sec. 5.4. 

As a test, we simulated the 3-state Potts AFM (using the same code, as explained in 
Sec. 4), so as to compare with exact results. (The 3-state model, in its representation as 
the BCSOS model, was also used for test simulations by Raghavan et al. (1997).) According 
to the theory in Sec. 2, limq_>o |qP(|^(q)P) = l/K. In Figure 5 we show |qp(|/i(q)p) vs. 
|qp for the 3-state simulation on a 64 x 64 lattice. Clearly this quantity docs approach a 
constant as |q| 0. Although the statistical uncertainty in each data point is ~ 10%, the 
uncertainty of the intercept in Figure 5 is only 1% according to a least squares fit. The 
fitted inverse stiffness constant is l/K = 1.92 ± 0.02, which is in excellent agreement with 
the known exact value of6/7r~1.91. 

5.1 Fourier mode fluctuations 

The expected behavior of the height components, in the limit — > and L — > oo, can be 
put in the form 

(|/^.(q)r)~^- (5.1) 
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When the height components are rough and described by a gradient-squared free energy, 
Eq. 3.3, then yi — and Q = (K"^) If we define height components to be 'smooth' either 
when (K^^),.,. = 0, or when the height fluctuations are bounded independent of system size, 
this is equivalent to j/^ > in (5.1). In such a case, the boundedness of height fluctuations, 
(|h — h'p) < const, implies, when inserted in (2.13), that operator correlations (O(x)O(x')) 
are bounded below as |x — x'| oo. One's usual picture of a 'smooth' phase is a background 
of long-range order, on which some height fluctuations with rapidly decaying correlations 
are superposed. That would imply g-independent behavior of (5.1) at long wavelengths, i.e. 
y = 2. 

Figure 6 shows log (g^(|/ii(q) p)) vs. log (g^) for the 4-state simulation on a 64 x 64 lattice. 
The slope of the graph at small q is yi/2 and its evident nonzero value implies that hi is 
smooth. A linear fit to the small-g data gave 

yi = 0.27 ±0.02, (5.2) 

where the uncertainty reflects the dependence on choosing an upper cut-off q for data points 
to be used in the fit. The data for h2 and (not shown) indicate similar values for y2 and 
^3, as expected from the symmetry among all even height components. 

As for the 'odd' height components, we find numerically for all q (7^ 0) that 

(|hodd(q)r) = 3(|heven(q)r). (5.3) 

This is seen in the nice overlap of plots of g^(|/ii(q)p) and g^(|/i4(q)P)/3 against (Figure 7). 
The corollary of (5.3) is that y^ has the same value for all height components. (The signif- 
icance of (5.3) is discussed in Sec. 5.2.) Hence all height components are 'smooth' and we 
expect the spins to exhibit long-range order. 

However, as noted above, the expectation for a smooth phase would be = 2, contrary 
to (5.2). Therefore, we carried out a finite size analysis to check whether the low-g behavior 
was a system-size effect. For nonzero q and finite L we expect (5.1) to generalize to 

{\k{q)\') = qy^-'Mq,L), (5.4) 
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where fi{q, L) is a scaling function that approaches Q as — > and L — > oo. Figure 8 shows 
a plot of g^(|/ii(q)p) vs. for 16 x 16, 32 x 32, and 64 x 64 lattices. The nice overlap of the 
three functions indicates that L in these runs is sufficiently large that fi{q,L) is a function 
of q only. 

The finite size of L also imposes a limit on the smallest nonzero wavevectors q which can 
be probed in a simulation. But the plot in Figure 6 appears rather linear (as already noted, 
its slope varies by only ±0.02 over a large range of q values), which would suggest that the q 
values accessible in the L = 64 simulation are small enough that fi{q, L) is close to /^(O, L), 
and hence that the value extracted for yi in (5.2) is close to the correct asymptotic one. 

5.2 Ratio of 'even' and 'odd' height fluctuations 

The ratio of 3 in Eq. (5.3) is nontrivial. It is certainly not true that |hodd(q)P = 3|heven(q)P 
for every allowed configuration; indeed. Figure 4 shows a case where |hodd(q)P is (essentially) 
zero while |hcvcn(q)P is unusually large for q = (0,±7r/L) or (±7r/L, 0). Nor could we find 
any hidden symmetry relating the 'even' and 'odd' height components. 

There is a simple explanation for the observation (5.3) if we postulate that 

(i) both even and odd height components are 'smooth', and 

(ii) the system is dominated by domains of a single ideal state, with small independent 
domains of the three other ideal states which are most similar to it. 

In the ideal-states graph, those neighboring ideal states have vertices surrounding that of 
the dominant state, forming 120° angles. (The same angles may be seen in both the even and 
odd projections of the ideal-states graph.) Consequently the height fiuctuations all occur in 
a rank-2 subspace of height space and in particular, hcvcn(q) is a linear multiple of hodd(q)- 
The linear ratio is in fact l/\/3 - the ratio of the even and odd components in Table 1 - and 
(5.3) follows immediately. 

In fact, if the fiuctuations occurred around a single ideal state, there would be cross- 
correlations between heven(q) and hodd(q)- However, thanks to the cluster update move, the 
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dominant domain is frequently being changed; thus cross-correlations are averaged to zero. 

In fact, we can even extend this argument to permit domains which are two steps away 
from 'AX' on the ideal-states graph. The second neighbor distances on the ideal-states 
graph are also in the ratio l/-\/3, so that the contributions of any single domain to heven(q) 
and hodd(q) ^-re still in the right ratio. However, the different domains do span a space 
of rank greater than 2, so the total hcvcn(ct) and hodd(q) are no longer linear multiples. 
Nevertheless, if the locations of the different domains are statistically independent, then 
cross-terms between different domains must cancel and one again obtains (5.3). Because 
small inter-domain correlations must exist, and there must be rare domains which are three 
steps away from the dominant one (third neighbor distances on the ideal-states graph are 
not in the right ratio), we conjecture that, in fact, there are tiny deviations from an exact 
ratio of 3 in (5.3). 

Thus, the observed ratio of 3 implies bounded fluctuations, which imply long-range or- 
der. Indeed, the only way we expect this behavior is when a finite fraction of the allowed 

configurations contain an AX domain that occupies a finite fraction of the lattice. 

5.3 System-size scaling of operators 

Say an operator 0(x) has correlations (O(x)O(x')) ~ r"'', and the corresponding order 
parameter is defined as Mq = ^ 0(x). Then simply expanding the double sum in {MqMq) 

X 

for a system oi N = sites gives 

(|MoH - L^-" - iV^-'?/^. (5.5) 

In particular, (|Mop) ~ N"^ in a system with long-range order. We measured expectation 
(5.5) for Mq —>■ M and Mq M^, corresponding to exponents t^m and rjs defined in Sec. 3.5. 

Figure 9 shows a log-log plot of (|M|^) from the 4-state simulations with three system 
sizes N. The data points fit to a line with slope 1.85 ± 0.01 (fitting error only!), i.e. 

(|M|2) ~ (5.6) 
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for N in the range simulated. 

We also measured the average total staggered magnetization for lattices of size 64 x 64 
and 120 x 120; connecting these two data points would imply 

{\Mf) ~ N^-^\ (5.7) 

Eqs. (5.6) and (5.7) respectively imply 

77m ~ 0.30 ± 0.02, 7;s~0.14. (5.8) 

5.4 Discussion: rough or smooth? 

The magnetization behavior reported in Sec. 5.3 is that expected for rough heights and 
critical spin correlations, which is in contradiction to our conclusion from Fourier fluctuations 
in Sec. 5.1. Which inference should we believe? We shall attempt to make sense of the data 
assuming a 'rough' state, and show it leads to several contradictions. 

We start with the h(q) data. The ratio of 3 in (5.3), which is our surest and most precise 
numerical result, implies that = ^K^. Then, in the notation of Sec. 3.5, it would follow 



that 



and hence 



VM^m^ rr""- (5.10) 



Roughness means that the curve in Figure 6 asymptotes to a constant as g — > 0; we may 
assume it does so monotonically, so the smallest visible value provides a lower bound on the 
stiffnesses, > 1.5, i.e. 

rj'^'' < 0.22. (5.11) 

The conclusion 77M = 77s in Eq. (5.10) is incompatible with the unequal values obtained 
in Eq. (5.8); furthermore, the value rju ~ 0.30 found from the magnetization data violates 
the inequality (5.11). 
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Most convincingly, since rf" jrp — 3/4 we find (from (3.15) and (3.18b)) that the locking 
exponent of Sec. 3.6 is r^iock = ^™ck ^ 0-2; far smaller than the critical value 4. Then (Jose 
et al. , 1977) the effective stiffness constant becomes larger and larger on longer length scales, 
so the system must lock into a 'smooth' phase. So, based on all the above arguments, we 
conclude that the system must be in a 'smooth' phase. 

We can check this picture by examining the system in real space. Long-range order, 
in which the system locks into a particular ideal state, means that a typical configuration 
should look hke an ideal state with some disorder: that is, there should be an excess of 
one particular spin label on one sublattice and a paucity of that spin label on the other 
sublattice, even in an infinite lattice. 

A measurement of this tendency can be gleaned from our results on the size distribution 
of XF-clusters (these were defined in Sec. 4.1.2). A substantial fraction of the encountered 
configurations had an Xy-cluster which included more than half of all sites; the frequency 
of such configurations decreases with lattice size, but is still quite appreciable at our largest 
lattice size, 200 x 200. Thus we cannot decide whether or not the frequency is zero in the 
thermodynamic limit. 

To further illustrate the tendency to long-range order. Figure 10 shows a typical config- 
uration (after equilibration) in the 4-state simulation. Note that roughly 2/3 of the sites in 
one sublattice have spin A whereas less than 3% of the sites in the other sublattice have spin 
A. In addition, all of the spin A sites in the ordered sublattice have the same height. The 
remaining sites are divided roughly equally among the remaining spin labels B, C, and D. 
These are the correlations expected in a long-range-ordered 'AX' type state. The fiuctua- 
tions from the 'AX' state (where the first sublattice has a spin other than A) appear to be 
domains of other ideal states that are one step from 'AX' on the ideal states graph; that is 
consistent with our explanation of the ratio of 3 in Sec 5.2. 

A final check of the long-range order hypothesis is to compute the total variance of each 
height component around its mean, which is {\zi{x)\'^) = ^q-^o(l^i(^)H- '^^^^ comes out 
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to be about 0.6 for the 'even' heights, which is also consistent with the picture just presented. 

As a tentative resolution of the contradictory features of our data, we suggest that our 
model may be close to, or even at, the unlocking (=roughening) transition. In that case the 
correlation length would be very large and it would be hopeless to detect the asymptotic 
behavior of the staggered magnetization in accessible system sizes. (We believe that the data 
in Figure 9 are not indicative of the asymptotic behavior as — > oo.) 

The anomalous power y — 0.27 extracted from the height fluctuation data might be 
explained if the proper functional form were 



This functional form is close to that for the 'rough' case, but strictly speaking the behavior 
would be 'smooth'. If we define an effective exponent ycs = 2— d log ((|h(q) p)) / d log q , then 
taking go to be the width of the Brillouin zone, we obtain ^eff — 1/ log {q/qo) ~ 1/ log(64) = 
0.24, which is comparable to the measured value. 




(5.12) 



32 



6 Conclusions 



We have found a spin model which has a height representation but still has long-range order 
in the ground state. It proved infeasible to detect the long-range order through studying 
fluctuations in the magnetization because the necessary system size is too large. However, 
by observing instead the fluctuations in the Fourier transformed height, we flnd conclusive 
evidence for long-range order through simulating systems of reasonable size. 

To discuss the possible flnite-temperature behavior of this model, and its relation to other 
models, let us display a Hamiltonian which has the CPAFM as its (degenerate) ground states: 

-H" = "2 XI ■^l(^a(x),a(x') - 2 XI ^25a(x),a(x") 
n.n. x,x' n.n.n. x,x" 

- X -'4(^fT(x),CT(x")<^CT(x'),CT(x"')) (6.1) 

plaq. x,x',x",x"' 

where the second sum is taken over second nearest neighbors, and the last sum is over 

plaquettes with x, x', x", x'" in cyclic order around the plaquctte. We must take the case 

J4 = — J2, making the plaquettes of type (ABAB) and (ABCB) degenerate, to make the 

CPAFM states all degenerate and justify their equal weights. 

At T > 0, typical conflgurations violate the CPAFM constraint at rare places; this 

creates dislocation defects in the height field; when h(x) is followed along a closed loop 

around that defect it is shifted by a height-space Burgers vector b (which is in the repeat 

lattice). The correlations of these defects are described by vortex- vortex (i.e. 'magnetic') 

type exponents in the Coulomb-gas picture (Nelson, 1983; Nienhuis, 1987). It can be shown 

that the specific heat and inverse correlation length behave as e"'^^"/^, where Eo is the 

energy cost of a defect, and the coefficient c is related to its Burgers vector and the stiffness 

constant(s) (Huse and Rutenberg, 1992; Read, 1992; Kondev and Henley, 1995a; Kondev 

and Henley, 1995b; Kondev and Henley, 1996).^ 

The best-known height models with critical ground states, in particular the 3-state Potts 
^This confirms that when Wang et al. (1989; 1990) fitted the specific heat of the 3-state Potts AFM to 
the form e^l^" , they should have obtained x = 1 (Ferreira and Sokal, 1995). 
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AFM, have stiffness constants K sucli tliat tlie defects are unbound, i.e., tlie system is in tlie 
same pliase as tfie two-dimensional Xy-ferromagnet above its Kosterlitz-Tliouless transition 
temperature. In such a case, T = is a critical point with the spin model being disordered 
at any finite temperature, since that means that the fugacity of the height defects, e~^°/^, 
is nonzero. A richer behavior is shown by the \axge-S triangular Ising antiferromagnet: the 
height field ensemble for the ground state is 'smooth' (i.e. the ground state has long-range 
order), and then as T" is increased there are two transitions at finite T, first to a critical 
state and then to a disordered one (Lipowski et al. , 1995). Since we find that the CPAFM 
model is 'smooth' at T = 0, it seems quite plausible that it undergoes analogous transitions 
at T > 0. 

The Hamiltonian (6.1) is equivalent to a Boltzmann weight (for each possible configura- 
tion of Potts spins) 

W ^ (1 - X5^(x),a(x')) 

n.n. x,x' 

X JJ^ [u-\-V (5o-(x),(t(x") + (^(t(x'),ct(x"')) + y(^o-(x),(T(x")<^CT(x'),(T(x"')] ■ (6-2) 

plaq. x,x',x",x"' 

The CPAFM constraints correspond to x = 1 and = in (6.2); the equal- weighted ensemble 
means also taking v — —y > Nienhuis (1990) studied the case where y = (see his Eq. 
(3)), which is a different weighting of the same constrained states. By mapping it to a loop 
model, he investigated this model with the number of Potts states as a continuous variable. 
It was found (see his Eq. (12)) that the model is critical when q < q^ — ?, (but would have 
long-range order when q > q^). Our Boltzmann weighting, with y < 0, disfavors the very 
flat 'ABAB' type plaquettes (relative to Nienhuis's weighting) and thus should be rougher, 
with qc> 3 in our case; this is all consistent with our conclusion that the model is smooth, 
but possibly is near the roughening transition. 

In Sec. 5, we speculated that the CPAFM model might be near the unlocking (roughening) 
transition. The 3-state Potts AFM on the Kagome lattice (Huse and Rutenberg, 1992) and 
the edge 4-coloring of the square lattice (Kondev and Henley, 1995a; Kondev and Henley, 
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1995b; Kondev and Henley, 1996) are known, in the ensemble where all configurations are 
weighted equally, to be exactly at their roughening transitions, i.e. the locking operator is 
marginal. In the latter case (|h(q)p), as measured in Monte Carlo simulations (Kondev and 
Henley, 1995b), shows a clean l/|qp behavior just like a non-marginal rough case. In the 
dimer-loop model (Raghavan et al. , 1997), one of the two height components shows clean 
l/|qp behavior indicating roughness, but the other component shows an intermediate power 
somewhat like our findings for the CPAFM model. 
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APPENDIX A: Dimensionality of height space 

In this Appendix we show that the height space has dimension 5 as asserted in Sec. 3.1. To 
do this, we first construct the most general mapping of a configuration to a height-hke vector, 
such that the height steps between two sites depends only on their respective spin values cr, 
a' and on the parity e of (say) the first site. This turns out to give 12 candidate height 
components. Then we show that only five of these are valid, nontrivial height components, 
and finally we show (by examples) that no further reduction of the dimensionality is possible. 

Note first that Az(cr, cr', e) = — Az((7', cr, — e) since the left- and right-hand sides are the 
height difference of the steps from x to x' and back again; thus, it suffices to define Az for 
e = -|-1. There are g(g — 1) = 12 possible combinations of (o. o') since neighbors cannot have 
the same spin value; thus all possible functions form a 12-dimensional linear space, labeled 
by (X, Y) e {A, B, C, D}^ with X y^Y. 

A candidate height component Cxy(x) is defined so Ad^xY = +1 if {<^, cr', e) = (AT, Y, +1), 
ACxy — —1 if (c", cr',e) — (F, X, — 1), and zero otherwise. In the CPAFM model, it can 
be checked that every component (xy is well-defined, by encircling each plaquette. (If we 
encounter an XY step, we must also encounter a YX step, e.g. in the plaquette XYXZ.) 

A.l Reduction to nontrivial height components 

Some linear combinations of the components Cxy(x) are bounded so the corresponding sub- 
spaces of the 12-dimensional candidate height space are trivial. Namely, let Po-,e(x) = +1 
when X is an even site and its spin takes the value a, and let p^^Ei^) = otherwise; let 
Po-,o(x) be the analogous function for odd sites. Then po._£;(x) = ~^a'^aCcTa'{^) (modulo 
a global additive constant), as can be checked by inspection of the possible /S.p^_E values; 
similarly po-,o(x) = + X^o-'/o- '^'^''^^"^^ (modulo a constant). There is one linear dependence 
among the 2g = 8 distinct p functions, since X^o.(Po-,£;(x) -l- Po-,o(x)) = 1. Thus there is 
^•^This assumption is justified since the height rule ought always to have the same lattice periodicity as an 
ideal state, which is a two-sublattice pattern for the CPAFM case. 
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a seven-dimensional trivial subspace of ( space; the remaining five dimensions, it can be 
checked, are spanned by the five height components defined in Sec. 3.1. (For the g-state 
CPAFM model, this construction would give — 3q + 1 dimensions.) 

A. 2 Necessity of all five components 

It remains to be shown that there is no further reduction of the height space dimensionality 
- that all five components are necessary. The easiest way to exclude the existence of a zero 
(or bounded and hence trivial) linear combination is by counterexample: we just exhibit five 
configurations, each of which has a net gradient only of the corresponding height component. 
We need only consider a one dimensional gradient. 

The bottom half of Table 2 shows the height steps corresponding to a four spin se- 
quence like {ABCD) (assumed followed by another A). Consider now the 10-spin sequence 
{ACBD){ABC){ACBD){ACB). If the first (ACBD) subsequence starts on an even site, 
then the second (ACBD) starts on an odd site; hence the even components of the steps are 
both (4, 0, 0)g but the odd components cancel. Also, the height step of (ABC) (beginning 
on an even site) is canceled (in every component) by the height step of (ACB) (beginning 
on an odd site). Thus the net height step is (8, 0, 0)^. A configuration, consisting of repeats 
of the 10-spin sequence, has a gradient of only the first even height component. By use of 
(ABCD) and {ACDB), we can construct a string of spins with a gradient of only the second 
or third even components. 

Furthermore, we could reverse the second copy of the 4-spin subsequence, obtaining 
(say) {ACBD){ABC){ADBC){ACB). This reverses the height step of the second 4-spin 
subsequence, so that now the even components cancel and the odd components add up. In 
this fashion we can construct a string which has a gradient of only the first or second odd 
height component. Thus, the five height components used in this paper are nontrivial. 

It should be noted that the above statement concerns only the allowed microstates of the 
system. It is possible - and indeed, likely in models with many height components - that 
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some of the components are "smooth". That is, although they can have arbitrarily large 
fluctuations (in an arbitrarily large system), the actual fluctuations are essentially bounded 
due to collective efi^ects. Thus, it may be that only a subset of the height components is 
described by a nontrivial elastic theory. 
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APPENDIX B: Proof for repeat lattice 

In this Appendix we rigorously show that two domains of the same ideal state always differ 
in height by a vector of the repeat lattice, from which it follows that the coarse-grained 
height of any ideal state is unique modulo the repeat lattice. This is done for both cases, 
the 3-state Potts AFM (see Sec. 2.2) and the 4-state CPAFM model (see Sec. 3.3). 

B.l Palindromic sequences 

In either case {q — 3 or q — A), we focus on the sequence of spins (7o(Ti(J2 . . . (J„_i(J„ connecting 
two sites. We do not care which path these spins lie along, since the height difference between 
the first and last site depends only on the sequence as in Eq. (2.6); or, for g = 4, on the 
sequence and on the parity eo of the initial site as in Eq. (3.5). We denote the height 
difference as Az(eo; (Joaia2 ■ ■ ■ cr„). Now 

Az(eo; (JocriC72 . . . aj^iajaj+i . . . On-xOn) 

= Az(eo; oqUxOi . . . crj--iaj) + Az{ej; (TjOj+i . . . Gn-iOn), (B.B.I) 

so we shall proceed by breaking all sequences down into a few kinds of subsequence. 

We call 'palindromic' a sequence of form a'^^'^a^^^ . . . a'^^^^^a^^^a'^^^^^ . . . a'^^^a^^^ (with j > 
1), i.e. the sequence has reflection symmetry about one spin. Of course, the corresponding 
spins both have the same parity. The net height change Az of any palindromic sequence is 
clearly zero, since (as can be checked from the height rule definitions) Az(e,;(7«(T(^+i)) = 
-Az( ej+i; cr'-^+^V*-*^). (In fact, this obviously must be true since the palindromic sequence 
also describes a path from the initial spin back to itself.) Therefore, using (B.B.I) the net 
Az of a sequence containing this palindrome is invariant if the palindrome is reduced to the 
single spin a^^^ . Note that such a reduction keeps invariant the parity of the length of the 
sequence. 
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B.2 Three-state Potts AFM case 

We are concerned with paths between sites in the same kind of domain - say, an AX domain. 
So we hmit ourselves to paths which connect sites of the same parity. Then the corresponding 
sequence of spins has an even number of steps, and it begins and ends with the same Potts 
state A. 

We can break each such sequence up into subsequences which begin and end in A; by 
(B.B.I) the net is the sum of those of the subsequences; furthermore, we can reduce 
any pahndromes within the subsequences without affecting this sum. Thus each term in the 
sum represents a pahndrome-free sequence beginning and ending with A and containing no 
interior A's; the only such sequences are A, ABCA, and ACBA, which have equal to 0, 
3, and —3, respectively. There must be an even number of three-step terms, since we chose 
the total sequence to have an even number of steps. Thus, finally, the net Az is a sum of an 
even number of ±3 terms, and hence is a multiple of 6, which completes the demonstration: 
the possible height differences A/i between two A B/C domains are precisely the integer 
multiples of 6. 

Analogous results hold for the macroscopic height difference between any two domains of 
(possibly) different ideal states. For example, Eq. (2.6) for A/i between an A B/C domain 
and a C/A B domain is simply that for some spin sequence of odd length that begins with A 
and ends with B. We can always decompose such a sequence into two subsequences, the first 
beginning and ending with A and the second beginning with A, ending with B and containing 
no interior A's. The first sequence has Az = (mod 6) and the second subsequence can 
be reduced to the pahndrome-free AB, for which Az = 1. Hence the macroscopic height 
difference between an A B/C domain and a C/A B domain is always congruent to 1 (mod 
6). 
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B.3 Four-state CPAFM case 

The method of proof for the 3-state model must be extended in this case since there are 
now infinitely many palindrome-free sequences beginning and ending with A and containing 
no interior A's. To address this, let us first replace every spin u*^'^^ which is not an A spin 
by the palindrome a^'^^Aa^''^; clearly this doesn't change Az. In this fashion the entire 
sequence is broken into short subsequences and we can write Az as the sum of terms of 
the form Az{e,AXA) or Az{e, AXY A), where X and Y are spins distinct from each other 
and from A. The palindromes AXA are produced whenever the original sequence had an 
A; such palindromes can be reduced leaving only the 3-step terms. Moreover, since they 
involve an odd number of steps, the parity must alternate between e — 1 and e = — 1 in 
successive terms. Thus, Az is the sum of terms of form Az(l, AXY A) + Az(— 1, AX'Y'A), 
where X and Y need not differ from X' and Y'. Constructing all possible terms of that 
form (using Table 2), one finds that Az must be a linear combination of the six vectors 
(4, 0, 0)e ± (0, 6)o, (0, 4, 0)e ± (-6, 0)o, and (0, 0, 4)^ ± (6, -6)o, which are the same as those 
in Eq. (3.6). (Conversely, it is also obvious that for any such linear combination, a sequence 
exists.) Clearly the Bravais lattice generated by hnear combinations of these vectors is the 
'repeat lattice' for the 4-state CPAFM. 
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APPENDIX C: Ergodicity of the update move 

Ergodicity requires that it be possible to apply successive update moves to any allowed 
configuration to obtain any other allowed configuration. Since the update move we used is 
reversible, it is sufficient to show that any starting configuration can be converted into the 
configuration with all A's on the even sublattice and all B's on the odd sublattice. 

Also, note that in an XF-cluster all the X's are on one sublattice and all the Y's are 
on the other sublattice. A corollary is that if a site with spin X and a site with spin Y are 
in the same sublattice, then the two sites are not in the same XF-cluster. Similarly, if two 
sites both have spin X but are in different sublattices, then they are in different XF-clusters 
for all F 7^ X. 

Now, consider an arbitrary initial configuration of spins. Pick any site and note its spin, 
X, and which sublattice it is in. If the chosen site is in the even sublattice and X ^ A, 
then interchange spins X ^ A within the XA-cluster containing the site; if the chosen site 
is in the odd sublattice and X ^ 'Q, then interchange spins X ^ 'Q within the XB-cluster 
containing the site; otherwise, do not perform an update move. In this way we set the spin 
on the chosen site to the value it takes in the 'target' configuration, which consists of A's on 
the even sublattice and B's on the odd sublattice. 

In the same manner, proceed one by one through the rest of the sites in the lattice, 
performing the appropriate update move, if necessary, to set the site's spin to the value it 
takes in the target configuration. The corollaries mentioned above ensure that once a site's 
spin has been set to its target value, it is not altered by the subsequent update moves used 
to set other spins to their target values; for, a site whose spin is set to its target value and 
a site whose spin is not set to its target value can never be in the same update cluster. For 
this reason, at the end of the process the target configuration has been reached. Therefore, 
the update move is ergodic. 

Note that this proof is valid whether or not periodic boundary conditions are used. In 
the case of periodic boundary conditions, however, it has been implicitly assumed that the 
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target configuration exists, which requires that the lattice have even numbers of rows and 
columns. As mentioned at the beginning of Sec. 4.1, this requirement is a natural one to 
include in any simulation since it is necessary for the ideal states of the model to exist. 

Also, note that this proof is equally valid for a standard (nonconstrained) Potts AFM so 
long as 'cluster' is appropriately redefined. In fact, the arguments here prove the ergodicity 
of the update move for any bipartite lattice since the target 'AB' configuration exists. This 
does not contradict the statement of Lubin and Sokal (1993) that the algorithm is not ergodic 
on 3m X 3n lattices when m and n are relatively prime because we require that m and n 
both be even. 
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List of Figures 

Fig. 1. (a) A ground state of the 3-state Potts AFM. On the left side is a portion 
of the A B/C ideal state; in the middle, a portion of the AB state (not an 
ideal state); on the right, a random state, (b) The corresponding height 
configuration. (Note that the heights are well-defined only after the height at 
one site is fixed.) 

Fig. 2. (a) The six ideal states of the 3-statc Potts AFM arranged in a circle, with 
heights as indicated, (b) The 'ideal-states' graph, with repeat distance 6, and 
operator components and m shown as a function of h. The operator P 
(not shown) would alternately take values +| and — i. 

Fig. 3. (a) Portion of a ground state of the 4-state CPAFM. The left portion shows 
the A B/C/D ideal state, while the right portion (separated by a vertical 
dashed line) shows a random state, (b) The corresponding configuration of 
even height component hi. (c) The corresponding configuration of odd height 
component h^. 

Fig. 4. (a) A 'roof starting configuration for the 3-state Potts AFM. The height 
variables in each quadrant have a (different) uniform tilt, which is the largest 
possible. This is also a valid 'roof configuration for the 4-state CPAFM model, 
(b) The corresponding configuration of even height component hi in the 4- 
state CPAFM. (c) The corresponding configuration of odd height component 
/i4 in the 4-state CPAFM. 

Fig. 5. Plot of |qP(|/i(q)P) vs. |qp for the 3-state simulation (size 64 x 64). The 
(nearly) horizontal line shows the best linear fit to the data. 

Fig. 6. Plot of |qp(|/ii(q)|2) vs. jqj^ (log-log) for the 4-state CPAFM model. 
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Fig. 7. This graph overlays plots of |qp(|/ii(q)p) and of |qp(|/i4(q)|^)/3, against 
|qp for the 4-state CPAFM model showing they cannot be distinguished. 
(Data are from a 64 x 64 lattice and, for clarity, only q || (±1, ±1) have been 
included.) 

Fig. 8. Scaling plot overlaying data on |qp(|/ii(q)p) vs. |qp for the 4-state CPAFM 
model, for various system sizes. (The largest q| values have been omitted.) 

Fig. 9. Mean square total magnetization (|Mp) vs. system size N (log-log plot). 
The line is the fit (|MH ~ N^-^^. 

Fig. 10. A 32 X 32 portion of a typical configuration found in simulation of the 
L = 64 system. 
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Spin pair 


Az 


AB 


(0,1,-1)e + 6(2,-1)o 


AC 


(l,0,l)E + e (-1,2)0 


AD 


(-l,-l,0)E + e(-l,-l)o 


BC 


(-l,l,0)E + e(-l,-l)o 


BD 


(1,0,-1)e + 6(-1,2)o 


CD 


(0,l,l)E + e(2,-l)o 



Table 1: 
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Spin sequence 


Az 


ABCA 


(-2,2,-2)E + e(2,2)o 


ABDA 


(2, 2,-2)^ + 6 (2, -4)o 


ACDA 


(2,2,2)E + e(-4,2)o 


ACBA 


(2,-2,2)^ + 6(2,2)0 


ADBA 


(-2,-2,2)E + e(2,-4)o 


ADCA 


(-2.-2.-2)^ + r(-1.2)o 


ACBDA 


(4,0,0)E + e(0,6)o 


ABCDA 


(0,4,0)E + e(6,0)o 


ACDBA 


(0,0,4)E + e (-6,6)0 


ADBCA 


(-4,0, 0)E + e (0,-6)0 


ADCBA 


(0,-4,0)^ + 6 (-6, 0)0 


ABDCA 


(0,0,-4)^ + 6 (6, -6)0 



Table 2: 



61 



